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Abstract 

We follow the emergence of quantum entanglement in a scattering event between two initially 
uncorrelated distinguishable quantum particles interacting via a delta potential. We calculate 
the time dependence of the Neumann entropy of the one-particle reduced density matrix. By 
using the exact propagator for the delta potential, we derive an approximate analytic formula for 
the asymptotic form of the two-particle wave function which is sufficiently accurate to account 
for the entanglement features of the system. 
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I. INTRODUCTION 



Quantum entanglement has become an important research topic in modern physics, 
not only because it exhibits the striking differences from classical concepts, but also since 
it is widely considered now as the fundamental resource in quantum information theory. 
Although the first famous paradox connected with entanglement ij] was presented in the 
context of observables with continuous values, the specific entanglement properties of such 
systems are less explored than those with discrete, e.g. spin states. 

Recent studies of continuous variable quantum systems focus on the emergence of 
bipartite entanglement in a scattering event of two interacting distinguishable quantum 



particles, which have no initial correlations 



Since this is a fundamental process in 



quantum physics, it is important to explore how does it gene rate quantum entanglement 
Some general features of this process were identified in [12| 



,Q. Refs. Is, 



Il0|,ll4 considered 

specific models on the scattering of ultracold atoms trapped in a harmonic potential 
well. Important results on the entanglement of colliding particles, modelled by Gaussian 
wave packets and interacting with different finite range potentials, were published in 



wave pa c 

flfly 



13 



15 



16|. 



In the present work we consider the explicit time dependence of entanglement in a 
quantum mechanical model of a collision process which itself creates the entanglement 
during the interaction between two particles that were independent in the beginning. To 
be specific, we assume a non-relativistic one dimensional motion with an attractive or re- 
pulsive delta potential between the particles. The evolution of the process is described by 
the explicit solution of the time dependent Schrodinger equation. Although this problem 
has been considered previously, here we derive analytic expressions for the post collision 
behaviour that incorporate the spread of the individual wave packets, by using the exact 



propagator for the delta potential 



13, 



18| . To quantify the entanglement we use the Neu- 



mann entropy [19|, and we present how it is built up during the collision. The asymptotic 
value of the entropy is then obtained from our analytic expression of the long time limit 
of the time dependent two-particle wave function. This study may find application e.g. 
in the experimental analysis of collision and recollision of atomic fragments following a 



laser induced dissociation 



20|. 
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II. INTERACTION OF TWO PARTICLES VIA A DELTA POTENTIAL 



The Hamiltonian of the system is written in terms of the position and momentum 
operators of the particles as 

H^2 = ^ + ^- Vo5 (Xi - X2) . (1) 
2mi 2m2 

For an attractive (repulsive) interaction we have here Vq > 0, (Vq < 0). We introduce, as 
usual for two-body problems, the operators: 



Xo = aiXi + a2X2, Po ^Pi + P2 (2) 
X =Xi-X2, P ^a^Pi-a^P^ (3) 

resulting in a sum of two independent Hamiltonians corresponding to the center of mass 
motion and the relative motion: 

H,2^Ho + H, ^0 = ^, H^^-Vo6{X). (4) 

Here m — mi +771-2 is the total mass of the particles, ai — rrii/m and // = mim2 /{fni +7712) 
is the reduced mass of the system. In the center of mass reference frame the expectation 
value of Pq; the total momentum of the particles, is zero, and the natural coordinate sys- 
tem is the one which has its origin in the expectation value of the center of mass operator, 
Xq. Then (Pq) = and (-^0) = for all times. We shall proceed by using coordinate wave 
functions and assume that initially the particles are described by a product of normalized 
Gaussians 

* (xi,X2,t = 0) = ^l{xi)^p2{x2) (5) 

(7-\/7r 

localized at distant points: around — oi = —a2a and 02 = ccio, as required by {Xq) — 0. 
Prom (Po) = we also have —k2 — ki —: q. In terms of the center of mass and relative 
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coordinates, Xq and x this wave function takes the form: 



{xi,X2,t = 0) = $ (Xo, X,t = 0) = ifc (Xo) (fir (x) , (7) 



where 



(^.(xo) = a-'/\-'/'e-^o/^'^\ (8) 
y^r{x) = (a/a)i/V-i/^e-°'("+'')'/2-^e'9" (9) 

are normahzed functions of Xq and x, respectively, a = ai + 02 is the mean value of the 
distance between the particles, and = aia2 = fi/m. 

The separability of the wave function in terms of the coordinates x and xn is due to the 

"nn — 



specific choice of the initial wave function p, ll2|, |l3|, |l5|, ll6|], where the variances of the 
positions of the individual particles (AXj)^ obey mi(AXi)^ = m2(AX2)^. In the more 
general case, one has a double sum of products of arbitrary basis functions in the new 
variables, which could still be transformed into a single sum in the Schmidt bases of the 
respective spaces (see Eq. flTB]) below). As implied by the linearity of the Schrodinger 
equation, the time evolution of the initial state could then be obtained by solving the 
problem for each term in the sum. 

The time evolution of the wave function in the coordinates Xq and x are determined 
by Hq and H independently, and they can be given by the respective propagators. For 
the free motion of the center of mass this is well known: 

Ki^o, yo, t) = (^-^ j exp [im{xo - yo)V2ht] (10) 
which yields the usual spreading Gaussian wave packet according to: 

$c (xo, t) = J K^{xo, Vo, t)(pciyo)dyo = 

= iV,exp[-xg/2a2] (11) 

where Nt = 7r~^/^(cr + iht/ma)~^^'^ and = o"^ + iht/m. The propagator for the delta 
potential Hamiltonian is more complicated, but still can be obtained _in a closed form. 
For the attractive case (Vq > 0) the propagator has been derived in 



while it turns 



out that the result is vahd for both signs of the potential and is given by: 



K{x, y, t) = Kl{x, y, t) + |e-T^e"'erfc {u) , (12) 



where 



/^K) n \x\ + \y\ / -f,./^ /ION 

The time dependence of the relative wave function can now be given as: 

<l>r{x,t) = J K{x,y,t)<^r{y)dy (14) 

It is not possible to determine $r {x, t) in a closed form given its initial value by 
we have to rely on numerical integration, but a very good approximate formula, valid for 
large times, will be given below. 

In order to consider the entanglement of the particles one makes the substitution 
corresponding to ([2]) and obtains a function of Xi and X2- 

^ {Xi,X2, t) = (Xo, t) (X, t) = ^daiXi + ^2X2, t)$r(Xl - X2, t). (15) 

which is not a separable state in the original coordinates xi and X2 any more. The above 
mentioned approximation for (x, t) will enable us to treat \E' {xi,X2, t) analytically, and 
consider explicitly the entanglement involved in it. 

III. REDUCED DENSITY OPERATOR AND ENTANGLEMENT 

In order to quantify entanglement in the state in (1151) one uses a measure that charac- 
terizes how much an actual two-particle wave function is different from a single product 
of two one-particle wave functions. In the context of quantum mechanics this was for- 



mulated first by J. Neumann [19|, based on the Schmidt decomposition [2l|] theorem. It 
states that for a square integrable function \E'(a;i, X2, t) of two variables, there exist a set 
of functions (j)k{xi,t) and ipk{x2-,i) which both form an orthonormal (but not necessarily 
complete) set in their respective Hilbert spaces, such that ^ can be written as a single 
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sum of their products: 



^'(Xi, X2, ^) = Vpk(t)(pk{xi, t)iJk{x2, t), (16) 

k 

The actual values of the p^-s are the simultaneous eigenvalues of the reduced density 
operators and ^2 describing either of the two subsystems defined by the hermitian 
kernel: 

gi{x[,xi,t) = J ^*{x[,X2,t)'^{xi,X2,t)dx2 (17) 



for system 1 and similarly for system 2. As it is shown in [19|, |2l| in more detail, ^1 
and §2 have a complete set of square integrable eigenfunctions forming a basis, 4>k{xi), 
ipk{x2) respectively, such that the corresponding eigenvalues pk are identical. The pk-s are 
nonnegative and form a discrete set, the sum of which is equal to unity. This also implies 
that the multiplicity of a positive eigenvalue must be finite. 

For quantifying entanglement it is natural to use the measure of randomness of the 
discrete probability distribution given by the pk-s in the Schmidt sum. Statistical physics 
tells us that this is best characterized by = —Y^pklnpk, which is just the Neumann 

k 

entropy belonging to the reduced density operator for each of the particles 



Sit) = -Tr [pi(t)ln/5i(t)] = -Tr [p2(t) lnp2(t)] (18) 

In order to calculate S one has to find the nonzero eigenvalues pk, that shall be time 
dependent during the collision. Fig. [T] shows numerical results about the time dependence 
of the quantum entanglement during the collision process, using atomic units. 



IV. APPROXIMATE PROPAGATOR AND ENTROPY FOR LONG TIMES 

The time evolution of the relative wave function, i.e. the second factor in f|T^ cannot 
be given in a closed form. We find an approximate analytical formula for the relative 
wave function and entropy using certain approximations for the propagator given by fll2p . 

The assumption that initially the particles are localized in a large distance from 
each other means that ipr{y) is different from zero only around —a <^ 0, therefore the 
contribution to the integral in f|T^ for y > can be neglected, and \y\ can be replaced by 
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FIG. 1: Time dependence of the quantum entanglement during the collision process. We plot 
the Neumann entropy (jlSp . computed numerically, vs. time measured in units of the classical 
collision time tc = ma/hq. The black solid curve is for rrii = m2 = 1, while the gray solid curve 
is for mi = 1, m2 = 10, as indicated in the figure. The dashed lines mark the asymptotic values 
of 5", given by formula (P5|) . We use atomic units, the parameters are: Vq = —5, a = 10, g = 
5, 0- = 1/2. 



—y in the propagator f|T2|) . Setting v = hq/m, which is the mean velocity of the relative 
wave packet, we can consider the asymptotic behaviour of the system for times much larger 
than tc = a/v, which is the time instant of the corresponding classical collision. We can use 
then the asymptotic approximation [22!: exp(M^)erfc(-u) = {y/n u)^'^{l — u^'^ /2 + 0{u~'^)), 
valid for large values of u. Keeping only the first term here we have then: 



K = K^ + i^f 



(47r/3)-i/2exp 



4/3 
4/3 



(19) 



The reduced (one-particle) density matrix in the form of Eq. fll7l) cannot be calculated 
from the propagator f|T9|) analytically. Therefore, we simplify it further by replacing the 
X and y variables of the propagator with the classical initial and final coordinate values, 



7 



y = —a and |x| = | —a + hqt/^\, in the pref actor of the exponential of Ki, but keeping the 
position dependence in the rapidly oscillating phase factor. This approximation is similar 
to the usual one in scattering theory, and yields the propagator 



K{x,y,t) = (47r/3)-i/2 



exp 



- yf 

4/3 



-1 



a + hqt/ /i| + a ^ 

ithg //i ^ 



4/3 



(20) 



Due to the presence of \x\ in the exponential of the second term, the integral (IT^ with 
the function i^riy) from ([9]) will split into two distinct Gaussians: one centered around 
X = —a + vt, this corresponds to the forward scattered wave, while the other around 
a — vt yielding the reflected wave. In this way we get the following asymptotic form of 
the relative propagator: 



K{x,y,t) 



2111 ht 

1 



c+ exp 



ifi{x - y f 
2ht 



+ c_ exp 



i/x(x + y) 
2ht 



21 



(21) 



where c_ = ";-^J/|7" ~ ^) ; c+ = 1 + c_. It is easy to check that for times larger than 
v/a these amplitudes coincide with the plane wave transmission and reflection coefficients 
for a delta potential with wave number q, which is the mean value of k in the initial 
relative state: 



lim c+ = T(g) := q/ (g — ig), lim c_ 



(22) 



The great advantage of the approximate form in ( 12T|1 is that it allows one to proceed 
entirely analytically and determine the total wave function, as well as the final value of 
the entanglement entropy in the system in a closed form. This is due to the emergence of 
Gaussian type integrals in f lT^ . which leads us to the approximate relative wave function 



^r{x, t) = T{q)(t)+{x, t) + R{q)(t)-{x, t) 



where 



(t)±{x,t) = Nt^/aexp 



a'^{±x + a — vty 



(23) 



(24) 



2(0-2 + m/m) 

The total two-particle wave function is now obtained with Eq. (fTSjl . and after some 
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algebra we obtain the result 



^{xi,X2,t) = [r(g)<l>i(xi,t)$2(x2,t) + i?(g)$i(-a;i,t)$2(-X2,t)] 



(25) 



with 




(26) 



and 



ht 



+ a^q[2{aj + ^ - ft/2] 



(27) 



ma- 



Here Oi = a2a, a2 = aia, the af/aj = (a^ + h^t^ / a'^m^) / aj describes the spreading of 
the respective wave packets, and Uj = hq/mj are the velocities of the particles in the 
corresponding classical problem. 

Fig. [2] shows $i(a;i,t) and ^i{—xi,t) in comparison with the two eigenstates of the 
computed one-particle density matrix, having the largest eigenvalues. 

The expression (1251) is the required approximation of the Schmidt decomposition, con- 
sisting of the two terms. The asymptotic value of the entropy of this entangled state is 
then 



The entanglement will be maximal for q = \g\, i.e. for center of mass momenta with —p2 = 
Pi = filVol/h, and has the value In 2. This explains why we have a larger entanglement 
for m2 = lOmi, when this condition is almost satisfied. 

V. CONCLUSIONS 

We have followed the emergence of the entanglement of two colliding particles that 
are independent initially and thus their wave function is a product state. In the case of 
the delta interaction potential and Gaussian initial states we have found an approximate 
analytic expression for the final entangled state, which is in good agreement with numerical 
results based on the exact propagator. Let us note that the wave function '^{xi,X2, t) is a 



S 



\nq)\ 



Hn\T{q)\^-\R{q)\^\n\R{q)\^ 




(28) 
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FIG. 2: Comparison of the asymptotic analytical results of the Schmidt decomposition, Eq. 
()25p . with the numerical method, for the case of (a) mi = m2 = 1, at t = 2tc, and (b) mi = 1, 
m2 = 10, at t = 2 tc- We plot the absolute value of the two eigenstates of the one-particle density 
matrix having the two largest eigenvalues, obtained numerically, with solid black (on the right 
side) and solid gray (on the left side) curves. The dashed black and dashed gray curves are the 
absolute values of ^i{xi,t) and ^i{—xi,t) functions, respectively, see Eq. We use atomic 

units, the parameters are: Vq = —5, a = 10, g = 5, a = 1/2. 
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kind of an EPR state [Ij in its original sense, i.e. in the coordinate space of two entangled 
particles. In [l| however the state considered is given by a highly singular delta function, 
while here \1/ {xi,X2,t) is square integrable during the whole process and has an especially 
simple asymptotic form. 
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